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q ■ Abstract 

O ' 
CM 

Assuming the S-matrix on noncommutative (NC) spacetime can still be developped 
perturbatively in terms of the time-ordered exponential of the interaction Lagrangian, we 
investigate the perturbation theory of NC field theory. We first work out with care some 
typical Green functions starting from the usual concepts of time-ordering and commu- 
^ ■ tation relations for free fields. The results are found to be very different from those in 

in 

^ ■ which can incorporate the new features of our results and which turns out to be the usual 

O ' 



the naive approach pursued in the literature. A simle framework then appears naturally 



time-ordered perturbation theory extended to the NC context. We provide the prescrip- 
tions for computing S-matrix elements and Green functions in this framework. We also 
D , emphasize that the naive seemingly covariant approach cannot be reproduced from the 

current one, in contrast to the field theory on ordinary spacetime. We attribute this to the 
phase-like nonlocal interaction intrinsic in NC field theory which modifies the analyticity 
c$ ■ properties of Green functions significantly. 
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1 Introduction 



The idea of describing coordinates in terms of noncommuting operators goes presumably 
back to Heisenberg and appeared already in a paper by Snyder |J in 1947, extended by 
Yang in the same year. It was associated with problems of treating hadrons which 
are extended objects and the assumption that there might exist a fundamental length. 
At that time renormalization theory was not yet well developped, but considered as a " 
distasteful procedure " 0. Even today this judgement might be shared by some people. 
In any case the opinion is widespread that at least in a theory which describes consistently 
all phenomena down to the Planck length the notion of spacetime undergoes some drastic 
change. The attempts of Connes and Lott |3[] aiming at a reformulation of the standard 
model in terms of noncommutative (NC) geometry started a new era which has been 
continued in the context of string theory Q. An analysis which brought the subject 
neatly into the context of somewhat more conventional quantum field theory has been 
provided by Doplicher et al 0. They showed that there exist representations of the 
coordinate operators such that for a class of minimal states in the state space of the 
system interactions can be written as the Moyal product of standard quantum fields as 
they have been introduced by Filk in another context |J. The road is then open to 
perturbation theory which has as a new element just this non-local interaction. 

Many calculations and model considerations have been performed in this vein and 
led to interesting new problems, most noticeably the ultraviolet-infrared mixing J7| and 
a potential violation of unitarity ||. They were all based on the assumption that up 
to modifying vertices by NC phase factors the Feynman rules are the usual ones, in 
particular that lines are represented by the conventional Feynman propagators. Doubts 
on this have been raised recently in a paper by Bahns et al. by reformulating in (p z 
theory the equation of motion in terms of the well-known Yang-Feldman equation JlC 



They arrived thereby at a solution which is manifestly Hermitian hence the theory must 
be unitary even if time-space NC is nonvanishing. This is, of course, in line with the 
general considerations of Ref. H, but does not yet explain why one arrives at violation of 
unitarity in the same model if one uses the Feynman rules of Filk ||, as shown in Ref. ||. 
Our study shows that the answer is very simple. Starting from standard free field theory, 
i.e. the standard commutation relations and multi-particle Fock space, one still formally 



defines the S-matrix as S = Texp 



. But when performing the contractions 
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according to Wick's theorem properly one can never combine the contraction functions of 
positive and negative frequency to the causal Feynman propagator. This arises because 
when time does not commute with space, the time-ordering procedure does not commute 
with the star multiplication either. The naive approach in terms of Feynman propagators 
is thus not well founded. 

In the next section we shall illustrate this statement by a detailed analysis of some 
examples. The picture on how to proceed in the general case will appear naturally in 
section 3. It turns out that the correct procedure of doing perturbation in NC field 
theory is just the time-ordered perturbation theory extended to the NC case. We shall 
spell out the prescriptions for computing S-matrix and amputated Green functions. We 
also provide an argument on how the seemingly covariant result in the naive approach 
cannot be reproduced from the time-ordered perturbation theory. We conclude in the last 
section. 



2 Analysis of some examples 

The basic quantity in quantum field theory is the Green function which determines the 
probability amplitude of S-matrix. We assume that perturbation theory for NC field 
theory can still be developped in terms of vacuum expectation values of time-ordered 
products of field operators so that the Green function is computed as a series expansion 
as usual, 

G(x u ■ ■ ■ ,x n ) = < 0|T ^i(xi) • • • (p n (x n )exp ij d 4 xC int ^ |0 >, (1) 

where £ int is the interaction Lagrangian. When one expands the above in £ int , T is meant 
to be taken before integration over spacetime is carried out. The time-ordering itself is 
also defined in the usual manner; for example, for free bosonic fields, it is 

T(ip 1 (x 1 )ip 2 (x 2 )) = t(x\ - xl)<f 1 (x 1 )ip 2 (x 2 ) + t(x° 2 - x1)ifi 2 (x 2 )(pi(x 1 ), (2) 

where r is the step function. (We reserve 8 for the NC parameter.) In this section we shall 
work out explicitly some Green functions in NC field theory using the above definitions 
and free field commutators. The differences to the naive approach of NC field theory will 
be clearly illustrated by them. 
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2.1 The vertex for the cubic scalar interaction 

Let us start with the following three-point function, 

G(xi,x 2 ,x 3 ) = J d A x 4 < 0\T ((p(x 1 )(p(x 2 )(p(x 3 )((f'k(p'k(p)(x 4 )) |0 >, (3) 
which would occur in ip 3 theory. The star product is defined as 



i 



exp(-^a»)/i(s)/ 2 (y) 



(4) 

y=x 



where x, y are the usual commutative coordinates and M „ is a real, antisymmetric, con- 
stant matrix characterizing the noncommutativity of spacetime, [x^, £„] = i6^ u . As in 
the usual field theory, the basic idea to work out the above quantity is to reexpress the 
time-ordered product in terms of the normal-ordered product (N) plus contraction terms 
arising from field commutators when interchanging field operators. Sandwiched between 
the vacuum state, only the completely contracted terms will contribute. In the usual field 
theory this is accomplished by Wick's theorem and a completely contracted term is given 
by a product of Feynman propagators. The same result is taken for granted in the naive 
approach of NC field theory. But as we shall show below this is actually not the case 
when time does not commute with space. 

In the example of eq. (|3|), there are 4! time orders but we only have to consider two 
groups of them, namely, x\ > x® > (x° and x®) and (x® and x° 4 ) > x® > x\, while others 
may be obtained by interchange of coordinate indices. For both of them, the following 
result will be useful, 

T((f 3 (f 4 -kif 4 -kip 4 ) 
= N(ip 3 (p 4 -kipi-k (p 4 ) 

+N (2D(0)<p 3 <p A + j d s ^(p 3 cp{x 4 - q) 



(5) 



+t u N (D u -k ip 4 -k ip A + ip 4 -k D u -k ip A + ip 4 -k ip 4 -k D 34 ) 

+ T 43 N (D43 * Lf 4 -k ip 4 + ip 4 -k D 43 -k Lf 4 + Lf 4 -k Lf 4 -k D 43 ) 

+T34 (2D(0)D 34 + J rf 3 /i q D(x 3 - x 4 + q) 
+r 43 (2D(0)D 43 + J d 3 fi cl D(x 4 -x 3 -q) 

where the star multiplication refers to x 4 . Some explanations are in order. We have 
denoted the arguments of functions by the indices of their coordinates when no confusion 
arises; e.g., (p 3 = ip(x 3 ), D 34 = D(x 3 — x 4 ) and r 43 = r(x° — x®). The above result is 
obtained by pushing annihilation (positive-frequency) operators to the right and creation 
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(negative-frequency) operators to the left using the basic commutation relation between 
them, 

(p(x) = ip+{ x ) + ip-(x), 
[ip + (x),ip-(y)} = D(x-y) ( 6 ) 

= J d 3 fj, p exp[-ip + ■ (x - y)], 

where rf 3 /i p = c? 3 p[(27r) 3 2_E , p ]~ 1 is the standard phase space measure with E p = Vp 2 + m' 2 , 
and = (\E P , p) (A = ±) is the on-shell momentum with positive or negative energy. 
Due to NC of the star product we have to take care of the order of ordinary functions in 
addition to the order of annihilation and creation operators when interchanging operators. 
For example, 

N(ip 4 -k if^-k ip 4 ) 
= [(f + * (f + * if + + if~ -k (f + * if + 
* ip~ * if + + ip~ * lfi~ * ip~ 

+(p~ -k ((p~±(p + ) + (<p~ * (p~)±(p + }(x 4 ) 
+2D(0)ip 4 + J d 3 ^<p(x 4 -q), 

where q^ = O^q 11 and * refers to the star product using —6^ u . The last but one term 
arises from contractions of the left and right fields with the middle one which is divergent 
as usual, while the last term comes from the contraction of the left with the right which 
is nonlocal because of the star product. These two terms are the origin of the third line 
and the last two lines in eq. (|^). These terms will finally contribute to the disconnected 
part of G(xi,X2,xs) and are thus ignored from now on. The second line in eq. (|^) does 
not contribute either when multiplied from the left by (px<f2 or from the right by (f2<pi 
and sandwiched between the vacuum. We are thus left with the fourth and fifth lines in 
eq. (|5p. Before we proceed, we would like to stress that for NC time-space (i.e., 8 0i 7^ 0) 
it is not permitted to change the order of time-ordering and the star multiplication since 
the former involves distribution functions of time and the latter contains derivatives with 
respect to time. In this general case, we cannot move T34 and T43 inside of the star product 
to form with D 34 and D 43 the Feynman propagator Dp(x 3 — x 4 ). This is the main source 
of difference in the current approach from the conventional one pursued in the literature. 
Only when time commutes with space, the two approaches become identical. 

Let us go back to the computation of G(x 1 ,x 2 ,x 3 ) and consider the time sequence 

Y 1^ <x 2 -""^ 3 45 

A = < 0\(piip 2 N (D34 * if 4 -k (p 4 + (p 4 -k D34 * <£>4 + If 4 -k <f 4 * D34) |0 > . (8) 
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Up to disconnected terms, only the combination (pfipZ ' ' ' <^4 ' " " ¥a ' ' ' contributes. Ap- 
plying repeatedly, 

¥>l"¥>2 ¥>I = Vlvtvt + D 24 ipf + ^14^2 , (9) 

we obtain, 

A = {D 14 ,D 24 ,D 34 }+, (10) 

where the star refers only to x 4 and we have introduced the completely symmetrized sum 
of the star products, 

{B±, B 2 , 5 3 }* = Ett 3 B n (i) * B n (2) * ^tt(3), (11) 

where 7r3 is the permutation of three objects. The above result is symmetric in x\ and x 2 
and thus applies as well to the time sequence x 2 > x\ > > x\. Including the case of 
x\ > a;", we have, up to disconnected terms, 

< 0\T ((pnp 2 (p 3 ((p -k (p-k(p)(x±)) |0 > , s 

= T 34 {£>i4, ^24, Aa}* + T43P14, Am, A3}*, for T 12 T 34 , 1 ' 

where Ti 2 T 34 denotes the time sequence {x\ and x°) > (x® and The opposite case of 
(#3 and x°) > (a;? and £ 2 ) is similarly computed to be, 

T34T12 : T 34 {Al, D 42 , -D34}* + T4 3 {Al, A2, As}*- (13) 

The other two pairs of cases are obtained by permutation of indices 1, 2, 3: 



^23^14 
^14^23 
^31^24 

T 24 T 31 



T 14:{D 2 4, -D34, -D14}* + T 4 i{D 24 , -D34, Al}*, 
Tl 4 {D 42 , A3, A4}* + T 41 {D 42 , D 43 , Al}*, 
^24^34, A4, A4}* + r 42 {-D34, D 14 , A2}*, 
7"24{ A3, Al, A4}* + T 42 { A3, Al, A2}*- 



(14) 



The above results can be combined into a compact expression. Since the connected 
contribution involves exclusively functions Dj 4 and D 4 j (j = 1,2,3), we hope that they 
are also accompanied exclusively by step functions Tj 4 and r 4 y This is indeed the case as 
in ordinary field theory. Actually all terms except {A4, D 24 , A4}* and {D 4 i, D 42 , A3}* 
are already in the desired form. There are three contributions to {A4, D 24 , A4}* arising 
respectively from the time sequences (x° and x 2 ) > x® > x 4 , (x 2 and x®) > x\ > x° A and 
(x 3 and x®) > x 2 > x 4 , whose union is identical to the sequence (x® and x 2 and x 3 ) > x\ 
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corresponding to the product of step functions T14T24T34. A similar combination occurs 
for {-D41, D42, -D43}*, so that 



< 0\T ((p^p^p * <p * <p)(x4)) |0 > 



(15) 



where Xj = ± defines the direction of time in the relevant function, 



'.ji 



T j4 , for Xj = + D x 3 
T±j, for Xj — — ' j4 



Dj4, for Xj = + 
D^j, for Xj = — 



(16) 



We would like to emphasize again that it is generally not permitted to move the step 
functions into the star products in eq. (0). Only when time commutes with space, we 
are allowed to do so and the result of the naive approach is reproduced, 



{Dp^xx - x 4 ), D F (x 2 - x 4 ), D F (x 3 - x 4 )}*, for 9 0i = 0, 



(17) 



where D F (xj — x 4 ) = t^D^ + t~ 4 DJ 4 is the Feynman propagator in coordinate space. 

To transform into momentum space we first work out explicitly the quantity in eq. 
(|T5|). Using the expressions, 



i j4 



iXj f°° A exp[-isj(x°j - x° 4 )] 



2tt 



dsj- 



Sj + ieXj 
J d 3 n Pj exp[-ip jXj ■ (xj - x 4 )], 



we have, for example 

Ai A 2 A3 ( f)Xi . f)\2 , f) 
'14 '24 r 34 ^14 * -^24 * u 
3 

= n 

d A pj 



3=1 
3 

n 

3=1 



a 3 
34 



Sj + ieXj 
iPXjipj) exp[-ipj ■ (xj - x 4 ) 



exp[-z(p 1Al ,p 2 A 2 ,P3A 3 



(19) 



(2tt) 4 "^ 

where we introduced p° = Sj + XjE p . and 

Px(k) 



expH(PUi,P2A a >P3As)]> 



A 



2£ k [A;°-A(£ k -2e)] 
= ^[AA; -^-^)]) ' 
(fci, /c2, ■ ■ ■ 5 k n ) = fcj A kj, 

i<3 



(20) 



with pAg = 1/2 d^p^q" '. The pair of parentheses introduced above has some properties to 
be used implicitly later on. It changes sign when the order of arguments is reversed. The 
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shift of sign in some arguments is identical to the shift of sign in the remaining arguments. 
It is important to note that the wedge product in the NC phase of eq. ( |19D involves only 
on-shell momenta of positive (Xj = +) or negative (Xj = — ) energy corresponding to 
propagation in the time direction of x® > x\ or x® < x\. Including all permutations of D 
factors in eq. (p~5[ ) which amounts to summing over permutations of the NC phase, and 
integrating over x 4 , we arrive at, 

G(x 1 ,x 2 ,x 3 ) 



3 

e n 

A1A2A3 j=i 



(21) 



(2n) 4 5 A ( Pl + p 2 + p 3 ) exp[-i(p 7r( i )A , (1) , Pn(2)K {2) , P*<?)K<& )]• 



x 

T3 



Transforming into momentum space is now straightforward, 



G(k u k 2 ,k 3 ) 
3 

G(x 1 ,x 2 ,x 3 ) 



= II J d 4x j exp(—ikj ■ Xj) 

3 (22) 
= (2tt) *5\h + k 2 + fc 3 ) E II /7> A 

A1A2A3 j=i 

X E eX PH(^(l)A. (1) , K(2)K {2V K(Z)K (3) )l 

where fcj's are the incoming momenta into the vertex. We have reversed the signs of 
variables Xj and A, used P-\(—k) = P\{k) and the property of the parentheses to remove 
the minus signs in the arguments of the NC phases. 

We make a few comments on the above result. First, the NC phases involve only on- 
shell momenta of positive and negative energies. For a given set of A^, the permutation 
sum of NC phases is 

2cos(k 1Xl ,k 2X2 ,k 3X3 ) + 2cos(k 2X2 ,k 3Xs , k 1Xl ) + 2cos(k 3Xa , k 1Xl ,k 2X2 ). (23) 

Note that the above does not simplify into 6 cos(ki Xl , k 2X , 2 ) etc. as it does in the naive 
approach where all kj Xj are replaced by kj. The reason is that while J2j hj — is always 
true this is generally not the case with kj Xj where kj is replaced by XjE^r. even if this sum 
vanishes for some configuration of Xj and kj, it cannot vanish for all configurations. For 
the case of identical fields considered here, there is even no such configuration at all due 
to kinematics. Furthermore, since the NC phases depend on time direction parameters 
Xj we cannot exhaust the sum over Xj by using iP + (k) + iP_(k) = iDp{k) with iDp{k) = 
i(k 2 — m 2 + ze) -1 being the Feynman propagator. These findings are completely different 
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1 A x fc o \ (24) 



from the naive approach. It is intriguing that such differences occur already at tree level 
in perturbation theory and we thus expect that the whole picture of perturbation theory 
will be altered. The differences arise from the fact that we are in general not allowed to 
interchange the order of the time-ordering procedure and the star multiplication. Only 
when 0Qi = 0, the star multiplication does not involve time derivatives and the NC phases 
in eq. fl22|) are independent of Xj, and then the differences disappear. 

The external lines in G may be amputated by multiplying by an inverse Feynman 
propagator (iDp)' 1 for each external line and noting that 

P x (k) = D F (k) Vx (k), 

2 V E^J 
We obtain the 1PI vertex for the above G, 

t(k 1 ,k 2 ,k 3 ) 

= (2,)V£^)Enh 3 (y] (25) 

X E eX PH(^(l)V)- k ^)K(2 V fc T(3)A w(3) )]- 

7T3 

2.2 The two by two scattering through cubic interactions 

To motivate our generalization in the next section, we consider the following four point 
function, 

G(x 1 ,x 2 ,x 3 ,x 4 ) 

= J d 4 x 5 J d 4 x 6 < 0|T(7Ti7r 2 X3X4(7r*o"*7r)5(x*o"*x)6) |0 >, ! " f> l 

which would arise from the Lagrangian of cubic interactions amongst real scalar fields, 

Ant = -^(vr *0-*7r)(x) - g x (x*0*x)(x)- (27) 

We have deliberately introduced nonidentical fields to avoid unnecessary complications 
due to many possible contractions amongst factors of identical fields, which just amounts 
to proper symmetrization of the NC phases as we saw in the above example. Our goal will 
be the S-matrix element of the two by two scattering tttt — > XX an d its crossed channels. 
It is clear that we should first contract the two a fields. We have, for x\ > Xg, 

< 0|T (7r 1 7T 2 X3X4(7r * a * tt) 5 (x * a * x) 6 ) |0 > 
= J d 3 fi p <0\---(7r*e- ip+ - X5 *7i) 5 ---(x*e +ip+Xe *x)6---\0>, ^ 



and for x® < x$, 

< 0\T (n^XsX^ * o" * 7r) 5 (x * o" * X)e) |0 > 



= J d 3 fi p < 0| • • • ( X * e- ip +- Xfl * x)e • • • (tt * e^'* 5 * tt) 5 • • • |0 >, 



(29) 



where the dots represent other fields appropriate to the time ordering and p refers to the 
a field, especially p° ± = ±E p = ±^/p 2 + mf. 

Next we consider contractions of 7r and x fields. Since [ni, Xj] = 0, the relative order of 
it and x fields is irrelevant; what is relevant is the order within the groups (x®, x®, x®) and 
(x®, x% Xq) respectively. Corresponding to x° > x® and x§ > x®, we have two possibilities, 
2I25T346 and T 346 Ti25. Let us study the 71 field contraction, 

< 0| • • • T(7r 1 7r 2 7r 5 * e ±ip +' X5 * tt 5 ) • • • |0 > . (30) 

There are 3! orders. For example, for x® > x% > the above becomes, up to disconnected 
terms, 

< 0| • ■ • 7^+7^ 7T 5 " * e ±ip +- X5 * 7T 5 _ • • • |0 > 

= < 0| • • • (D 15 ir£ + D 25 irt) * e ±ip +- X5 k 7r 5 " ■ ■ ■ |0 > (31) 
= (£>is * e ±ip +^ ★ D 25 + (1 <-> 2)) < 0| • • • |0 >, 

where -D15 and -D25 refer to the 7r field and * refers to x$. The above is symmetric in 

x\ and x 2 and thus actually corresponds to the time order specified by T15T25. The other 

time orders can be similarly computed. Their sum gives the complete result for all orders, 

< 0| • • • T(7ri7r 2 7r 5 * e ±ip - X5 * tt 5 ) • • • |0 > 
= E^ 1 4 2 Pi5 1 *e ± ^^ J D 2 A | + (1^2))<0|---|0>. (32) 

Ai A2 

The last factor in the above is precisely the one for the x field contraction and is similarly 
computed. We thus have, 



< 0|T (71-^2X3X4(71- -ka-k 7r) 5 (x -ka-k x)&) |0 > 



dVfc £ ^^^^ 



X 



* e-^+^s * + (1 <_> 2) £>*g * e +ip+ - X6 * + (3^4) 
+ ( same as above except r 56 — > r 65 , x 5 6 — > —2:5,6) 

C ^ 3 / i P £ E r 15 T 25 r 36 r 46 4 T 56 



(33) 



/ 



x 



{A,} A 



D^l * e"^ 5 * L^ 2 + (1 <-> 2) L^ 3 * e +ipx - X6 * Z^ 6 4 + (3^4) 



where -DgJ and D^g refer to the x field and the star in the second factor is with respect 
to Xq. In the second equality we have made the shift p — > — p for A = — . 



10 



Using the representations as shown in eq. fllTf ) and the same trick that led to eq. (|l9l) , 
we make the above Green function ready for transformation into momentum space. For 
example, 



4 

n 

3=1 

xe 



Ai A2 A3 A 4 A 
I ' 
4 , 



r 15 r 25 T 36 T 46 T 56 



d Pj 



(2k) 



d 4 p 



iP\(p) 



(2kY 



(34) 



x exp[-i(p 1Al , -p A ,p 2A2 )] exp[-z(p 3A3 , +Pa,P4A 4 )] 



We can now integrate over £5 and xq, which results in two factors of 6 functions, then 
transform into momentum space, and sum over all terms, 



G(h, k 2 , k 3 , fc 4 ) 
4 



n 

3=1 



d x^e 



G(xi,x 2 ,x 3 ,x 4 ) 



(35) 



= [2k)h\y. h) e e n K (%)] ^aCp) 

i {A,} A j 

x (exp[-z(A; 1Al , -p A , fc 2A2 )] + (1 <-> 2)) (exp[-i(k 3Xs , +p x , k AXi )} + (3 «-> 4)) 
= (2vr)V(E^) E Ell [^(%)] ^Pa(p) 

i {Aj} A j 

x2cos(fci Al , -p A , A;2A 2 ) 2cos(A;3A 3 , +Pa, &4A 4 ), 
with p = fci + k 2 = — k 3 — k 4 . Be careful that fci Al + k 2 \ 2 ^ P\ ^ — k 3 \ 3 — ^4A 4 - For 
comparison, in the naive approach the above would be 



(2n)H 4 (j:iki)U j [iD F (k j )} i D F (p) 2cos(h,k 2 ) 2cos(k 3 ,h) 



(36) 



which according to our preceding analysis is correct only for 9^ = 0. The 1PI function is 
obtained by amputation, 

t(k 1 ,k 2 ,k 3 ,k 4 ) 



(27r)V(E^)EEIlK(^)] ^(P) 

i {Aj} A j 

x 2 cos ( fci Ai , -Pa , &2A 2 ) 2 cos ( k 3Xa , +Pa , &4A 4 ) • 



(37) 



Let us now extract the S-matrix element for on-shell particles from the above 1PI 
function. We take the example of tttt — > XX scattering. This means, k\ = +Ey ll , k® = 
+E\f_ 2 , k® = —Ek 3 , k® = — -Ek 4 - Thus only one term in the sum over {Xj} contributes due 
to r)-(ki) = etc. Including the coupling factors as well, the transition amplitude is 

iT (n(k 1 ) + n(k 2 )^ X (h) + x(h)) 
= (2vr) 4 5 4 (A; 1 + k 2 - k 3 ~ h) E iPx(p)(-ig*)(-ig x ) 



x2cos(k 1+ , -p A , k 2+ ) 2 cos(-k 3+ , p A , -k A+ ) 
= [2K)H\k 1 + k 2 -k 3 - h)iA, 



(38) 
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where p = k\ + k 2 = k 3 + k 4 and A is the amplitude with the usual normalization as 
computed from Feynman diagrams in ordinary quantum field theory. Note also that we 
have reversed the signs of k 3 and k 4 so that k 3 ^ — > — k 3+ , k 4 _ — > — k 4+ . More explicitly 

T(n(k 1 )+n(k 2 )^x(h)+x(h)) 

= -g^gxi^Y^iki + k 2 -k 3 - k 4 ) 

2 cos(/ci+, -p+, k 2+ )2 cos(-/c 3+ , p+, -k 4+ ) 

— (dyj 



X 



+ 



2£ kl+kl (hi + k° 2 - E kl+kl + ie) 
2 cos(-/c 3+ ,p_, -k 4+ )2 cos(k 1+ , k 2+ ) 



2E k3+k4 (-k° 3 - kl - E k3+U + ie) 

We make a few remarks concerning the S-matrix calculation. First, the crossed chan- 
nels of the above process may be obtained similarly. For example, for nx — *• ^X scattering, 
we may choose k\ = +E kl ,k 2 = —E k2: k^ = +_E k3 ,A; 4 = —E k4 . Second, from the above 
detailed analysis it is clear how to calculate the most complicated case of identical par- 
ticle scattering through their self-interactions. We should include all possible Feynman 
diagrams and for each of them employ the same analysis which just amounts to more 
symmetrization at the vertices with respect to identical fields. In this way we get the 
following contributions to the amputated four-point Green function of the tt field at the 
lowest level in £ int = — g n Tr * tt * tt, 



t(k 1 ,k 2 ,k 3 ,k 4 ) 



4 



= -<£(2tt) 4 <J 4 (5» EEII K-(%)] (A s + A t + A u ), 

i A j=l 



(40) 



where A s ,A t ,A u are from s— , t—, u— channels respectively, 

A s = iPx{pa) X] ex PH(^iAi, hx 2 , -Ps\)] X]exp[-i(/c 3A 3,/c 4A4 ,+p sA )], 

7T3 7T3 

A t = iPx(pt)^exp[-i(k 1 x 1 ,k 3 x 3 ,-p tx )]^ex.p[-i(k 2 x 2 ,k 4 x 4 ,+Pt\)], (41) 

7T3 7T3 

= ^A(p«)Xl eX P[^( /i; lAi ) /i; 4A4,-P«A)]5I eX P[ - ^( /i; 3A3,^2A 2 ,+P«A)], 



7T3 T3 

2 exp[-i(fci Al , fcsAs, -Pu)] X! 

7T3 7T3 

^ expf-^fciAu A; 4A4 , -Pu\)\ X 

7T3 7T3 

with p s = k\ + k 2 ,p t = ki + k 3 ,p u = k\ + fc 4 . Here 7r 3 refers to the 3! permutations of the 



three momenta appearing in each factor of the NC phase sums. It is straightforward to 
project the S-matrix element from the above which we shall not write down. And there 
is also no problem to extend to more complicated interactions like ip A . Finally, if our aim 
is restricted to the S-matrix for on-shell particles, we may proceed more directly from 
the expectation values of the S-operator. For the above example, we need compute the 
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following quantity, 

< XX\S\im > 

= {-ig*)(-i9 x ) J dAx 5 J dAx 6 ( 42 ) 
+ higher orders, 

which just corresponds to a special assignment of the time order in eq. (p6j ), namely 
7i"i and 7T2 in the far past, X3 an d Xi in the far future, and others in between. This is 
precisely the contributing part in the above analysis. It is thus no doubt that the results 
for the S-matrix coincide. The main advantage of coping with the Green function for this 
purpose is that we may project all physical processes from the same Green function. 



3 Generalization 

The structure shown in eq. ([39|) looks familiar to us and is very suggestive. Actually it 
is nothing but the " old-fashioned ", time-ordered perturbation theory UTTfl ||12|| modified 
properly to NC field theory. This fits also on a somewhat more formal level to the 
pragmatic point of view which we assume here. The S-operator maps on the one hand any 
prepared incoming state onto the respective outgoing one, but this can also be considered 
as the time tranport of this incoming state into the outgoing one. As long as we can 
represent it as a time-ordered exponential with a Hermitian exponent times the pure 
imaginary unit % we have formal unitarity and thus satisfy the first requirement for true 
unitarity which means conservation of the transition probability in the sense of quantum 
mechanics - to which we shall return in a separate publication. 

The time-ordered Feynman diagrams for the above example are depicted in Fig. 1. 
In the language of the time-ordered perturbation theory, a physical process is virtualized 
as a series of transitions between physical intermediate states that are sequential in time. 
The transition amplitude is weighted by the interaction vertices which are evaluated for 
on-shell momenta if they depend on them and by the energy deficit of the intermediate 
states. Realizing this, it becomes obvious how to proceed in the general case. For further 
applications, we give below the prescriptions for computing the on-shell transition matrix 
T at some fixed order in perturbation, which are readily generalized from the ordinary 
ones [11 [12 . The additional piece for general off-shell amputated Green functions will be 



described later on. 
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(a) 



(b) 



Fig. 1: Time-ordered diagrams corresponding to eq. fl39|). Time flows upwards. 

1. Draw all Feynman diagrams for the process under consideration. For each Feynman 
diagram draw all of its time-ordered diagrams. Only the time order of interacting 
vertices is relevant and indistinguishable time orders are counted only once. Each 
time-ordered diagram is computed by putting together the following factors. 

2. Associate with each internal line (with a spatial momentum k) a phase space integral 



3. Associate with each vertex v, which is formed by internal lines j and external ones 
e of incoming spatial momenta q a and which has the interaction pattern in the 
Lagrangian —gipi * ip2 * • • • ipm an interaction factor g exp[— i{q\\ x , q2\ 2 , • • • , QnX n )] ■ 
X a = +(— ) if the vertex v is the later (earlier) end of the line a. The initial (final) 
particles e are always counted as earlier (later) than the vertex v. Symmetrize the 
above factor with respect to identical fields. Impose spatial momentum conservation 
at the vertex v by multiplying (27r) 3 5 3 (X)a Qa)- 

4. Associate with each intermediate state occurring between two sequential vertices 
(earlier v\ and later t> 2 ) a factor of energy deficit, [J^ei^Pe) ~ J2jE(kj) + ze] _1 . 
Here Y^ei^Pe) * s the algebraic sum of the zero-th components of external momenta 
entering (+) or leaving (— ) the diagram before and including the earlier vertex v±. 
E(kj) is the on-shell positive energy of the j-th line contained in this intermediate 
state. 

5. Multiply by a global factor of — 2Tr5(J2 e {^ z Pe)) an d a symmetry factor 1/S which 
excludes that of indistinguishable diagrams mentioned above. 

The above prescriptions would be precisely the same as obtained in ordinary field 
theory if we could interprete the vertex factor as a kind of numerator arising from spin. 
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In ordinary relativistic field theory we can recast the time-ordered perturbation theory into 
a covariant form in terms of Feynman diagrams. So, it is tempting to ask why this is not 
possible in NC field theory. Of course, Lorentz invariance is lost at the very beginning and 
it is not guaranteed that a seemingly covariant formalism exists and is equivalent to the 
time-ordered one if it does. But this is not the whole point. As far as Feynman diagrams 
are concerned, we can always treat 9^ as if it were a Lorentz tensor and there will be 
no problem if we do not use any special reference frame for calculation since we could 
not return back by a transformation afterwards |13|]. We could also consider 6^ v as some 
background field and assign to it a transformation law so that the above consequence still 
applies. In the following we present an argument that in NC field theory formulated via 
time-ordered perturbation theory one cannot reproduce the seemingly covariant results of 
the naive approach. Our time-ordered version seems however to be a safe starting point 
as far as quantum mechanics still applies to NC spacetime. A key element of it is the 
highly nonlocal character of NC interactions. 

Let us first recall briefly how to shift from the covariant perturbation theory to the 
time-ordered one in ordinary field theory. For a detailed account of the topic we refer the 
interested reader to Ref. ]12] for a nice presentation. For this purpose, we consider the 
contribution from a Feynman diagram to a general, unamputated and connected Green 
function in momentum space. One first expresses the 5 function of the zero-th components 
of momenta at each vertex in terms of a time integral. Then, one rearranges the product 
of time integrals thus obtained in a time-ordered way. This is followed by integrating over 
the zero-th momentum components, which is typically of the following form, 



Here t is the difference of the time variables introduced above between the two vertices 
connected by the line carrying momentum in the same direction of t. f(ko) is usually 
a polynomial of finite degree and analytic in the complex ko plane. The above integral is 
evaluated using contours. For t > (<)0, one closes the contour in the lower (upper) half 
plane picking up the residue at k = E k — ie (k = —E^ + ie) with the result, 



A crucial condition for the above manipulation is that f(ko) must not blow up at the 




(43) 




(44) 
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lower (upper) infinite semicirle in the complex plane faster than e~ tkot decays. Finally, 
one completes the time integrals sequentially and arrives at the result in the time-ordered 
perturbation theory. Now let us try to do the opposite in NC theory from the time- 
ordered perturbation theory to the covariant one by turning around the above procedure. 
In this case the function f(±E k ) is an NC phase which is essentially an exponential 
(superposition) of the form exp(±ii? k A; /0 ), where k'° = 0i k' z with k' being a spatial 
momentum of some other internal or external line. If the above procedure were reversible, 
the corresponding function in the complex A; plane would be something like exp(ikok'°) so 
that the naive result might have a chance to be recovered. But this is impossible because 
it is not guaranteed that it increases slower than e~ lkot decays. Actually whether it decays 
or blows up depends on the sign of k'° which itself changes with k'. This thus interferes 
with the above contour integration. One may argue that we may shift t to absorb k'°. 
This is again not legitimate since t is a time difference and doing so simply spoils the 
time-ordering procedure which is a key bridge to relate the two formalisms. Furthermore, 
a connected diagram has certainly more than one line; the above shift of t's, if it worked 
at all for one of them, would also interfere with each other making the trick totally useless. 
The above argument fails only if at tree level k'° is an external momentum and happens 
to vanish. But this is a very special kinematic configuraton if possible at all. We surely 
cannot rely on this in favour of the naive approach. From this analysis it is also clear that 
the main obstacle originates from the nonlocal exponential interaction that is intrinsic in 
NC field theory. 

Finally we extend the above prescriptions to amputated and connected Green functions 
by adding the following rule concerning external lines for time-ordered diagrams. It should 
also be applied to each individual diagram in which the connection of external lines to 
vertices is fixed. 

6. Multiply by a factor of rj\ e (p e ) for each external line with incoming momentum p e 
and time direction parameter A e which is +(— ) if it connects to a(n) later (earlier) 
vertex. This same A e also appears in the preceding vertex factors (where it takes 
one of the signs for S-matrix). Sum over the set {A e }. 
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4 Discussion and conclusion 



Based on the assumption that a time-ordered expansion of a formally unitary time evolu- 
tion operator is a good starting point also on NC spacetime, we studied perturbative NC 
field theory which turns out to be the time-ordered perturbation theory adapted properly 
to the NC case. This was achieved by a detailed analysis of some exemplifying Green 
functions which we worked out with care, and then extended to the general case. We 
found no obstacles in implementing noncommutativity in perturbation theory whether 
time commutes with space or not. We provided prescriptions for computing S-matrix 
elements and amputated Green functions. 

Although we only treated scalar fields in this paper, we expect no problems with 
spinor fields as it is the case on ordinary spacetime. Since the spinor effect amounts to 
an additional numerator associated with a propagator, it plays a similar role as a vertex, 
namely, the momentum contained therein will be on-shell with positive or negative energy. 
The situation is more complicated for gauge bosons. But again as in the usual theory 
there should be no problem at least in the 't Hooft-Feynman gauge. Our method also 
applies to any dimensions. 

The NC perturbation theory thus obtained is already different at tree level from the 
naive approach followed in the literature. The interaction vertices involve only on-shell 
momenta of positive or negative energy of participating particles. The basic quantity 
connecting vertices is not the causal Feynman propagator but the individual propagation 
functions of positive and nagative frequency. These elements are naturally incorporated 
in the framework of the time-ordered perturbation theory. In contrast to the ordinary 
field theory, it seems impossible to recast the NC time-ordered theory into a covariant 
form as has been assumed in the naive approach. We attributed this difference to the 
highly nonlocal character of phase-type NC interactions which has significant impact on 
the analyticity properties of Green functions in the complex energy plane. 

Since the whole picture for perturbation theory has been changed, we expect some of 
the phenomena found previously in the naive approach will also be altered. Amongst them 
we would like to mention briefly the unitarity issue which we will detail soon in a separate 
paper. Since the NC phase now involves only on-shell momentum, it is independent of the 
zero-th component of a generally off-shell four- momentum. The analyticity properties of 
Green functions in the complex plane of the zero-th component will thus be very different 
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from that in the naive approach. Furthermore, the right-hand side of the unitarity relation 
for Green functions is also modified due to the change of vertices, which does not seem 
to have been noticed thus far. Considering all of this it is quite reasonable to expect 
that unitarity will be practically preserved as it is formally built in the time-ordered 
perturbation theory. This may have some hints for the interplay of string and field theory. 
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